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Abstract 

I estimate maximum thermal conductivity k of a perfectly aligned bundle of single wall carbon 
nanotubes. Each row of aligned nanotubes has a discrete structure. It consists of segments of 
nanotubes with length L. The spacing between the segments block the phonon path through the 
row. Only the scattering due to the finite length of the segments is taken into account. The result 
is that the 'effective" mean free path is of the order of L/7. For 1 micron tubes (10,10) we get 
maximum value of k ~ 300W/m K at room temperature. This result is in a reasonable agreement 
with the experiment by Hone et al. assuming that in their samples L ~ l^um 
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INTRODUCTION 



It is expected that a single wall carbon nanotube(SWCN)is a very promising object for 
creation of metamaterials with a high thermal conductivity (TC) 0| • The first reason for 
this expectation is that the carbon-based materials, like diamond, have the largest known 
TC and the second reason is a molecular perfection ^ of the (SWCN's). However, to the best 
of my knowledge, the highest TC ever observed in SWCN's bundles at room temperature 
is about 220 W/mK and it is ten times smaller than TC of the natural diamond j3|. This 
highest result has been reported by Hone et al. Q] for a bulk sample of magnetically aligned 
nanotubes. The aligned SWCN's form a bundle where all tubes have a preferable orientation 

n 

in some direction. Hone et a/. [4] show that the TC of the aligned SWCN's is strongly 
anisotropic with the largest value in the direction of the alignment. The enhancement of the 
TC due to the alignment has been observed also by Zhou et a/.j5| and by Choi et a/.|6|], but 
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Q are 



the absolute values of the reported TC have been significantly smaller than in Ref. 4 
are many theoretical works on TC of the SWCN's. Some computational ones0, 0, 2 
made by molecular dynamics simulations. The results of these simulations have different 
values and different T-dependences. They predict mostly very high values of the TC. We 
think that the main problem of all these works is a small size of an array that can be 
simulated. There are also some different analytical approaches to the problem jlll. Il^ and 
beautiful reviews 3, 14. isj]. 

The purpose of my work is to estimate the maximum TC value of aligned nanotubes 
taking into account that they consist of segments with a finite length. It is well known, that 
tubes in ropes are not infinitely long, but have brakes, because each method of synthesis is 
able to create separated tubes of only a certain length. It is believed that this length is of 
the order of a few /im (See Ref. [16] and references therein) Then due to Van der Waals forces 
the tubes stick together and create bundles where the end of tube has no chances to make 
a strong chemical bond to the end of neighboring tube. There are many experiments that 
show, that tubes inside bundles have free ends. The idea of my work is to argue that this 
effect may be responsible for the relatively low TC as compared to the crystalline carbon 
materials. 

I consider a bundle of nanotubes perfectly aligned in x direction. Each segment of a 
nanotube has a finite length with an average value L. The nanotubes are organized in an ideal 
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triangular lattice with 6 nearest neighbors 



17l |. The cross section in a plane perpendicular to 



the nanotubes is shown in Fig. ^ (a). The cuts in each line of the nanotubes have random 
positions. Thus, on the length of each segment there are in average six cuts of its nearest 
neighbors. A homogeneous interaction between infinite tubes does not cause the loss of the 
phonon momentum. However, a phonon flux has to overcome the openings between the 
segments at the termination points of each nanotube segment. I assume that this openings 
are so large that a jump of a flux occurs with an assistance of all six neighboring rows of 
the tubes as it shown in Fig. ^ (b). Slightly different mechanism of momentum scattering 
appears in a given nanotube ("0") if one of the neighboring nanotubes has a termination 
point as shown in Fig. ^c). 

Having in mind to get a maximum estimate of the TC, the propagation of heat flux 
Q between the scattering points I assume to be ballistic. Quick phonon exchange at the 
scattering points leads to a thermalization of symmetrical parts of distribution functions of 
phonons with temperatures that are determined by values of effective thermal resistances 
between scattering points. It should be noted that my calculations cannot be applied directly 
to multiwall carbon nanotubes. 



SCATTERING PROBLEMS 



For simplicity an elastic scattering only is considered. This means that phonons generated 
in neighboring nanotubes not only have the same frequency, but belong to the same mode 
as the incident wave. It is not necessarily true but it should not lead to a big mistake 
due to small matrix element of the inelastic scattering. On the other hand, without this 
approximation the calculation of the TC would be a mess because one cannot consider each 
mode independently. I assume a linear dispersion law for the modes under study and argue 
that it is also not important if the interaction is large. 

Consider a propagation of some mode along the nanotube "0" as shown in Fig. (a). 
There are two basic scattering geometries for this mode. One of them appears at the 
terminational points of this very nanotube. It is shown in Fig. H(b). Only two neighboring 
nanotubes are shown, however all six are participating in the scattering. The interaction 
between the neighboring nanotubes is shown by springs in Fig. ^ The energy of this 
interaction is expected to be larger than Van der Waals interaction because divP at the 
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FIG. 1: (a) The cross-section of the bundle that show nanotube "0" and its nearest neighbors, 
(b) The first scattering problem. The cross-section by the plane of nanotubes 2-0-5. The wave 
incident from 0_ reflects backward and transmits through the opening into O4. with simultaneous 
excitation of the waves in all six neighboring tubes, (c) The second scattering problem. The wave 
incident from 0_ and is scattered by the cut in tube 2. It reflects backward, transmits into O-j. and 
excites waves in tube 2+. 

termination point might be non-zero, where P is a polarization in the wave. The second 
scattering geometry (Fig. ^c)) appears if one of the nearest neighbors of nanotube has a 
terminational point. In both cases only an elastic interaction is considered. 

In the first scattering problem the acoustic wave in the nanotube is described by a 
following equation 



(1) 



(2) 



where s is velocity of the mode, is a characteristic of the springs, i = 1 . . . 6, f/=p is the 
displacement of a tube "0" at x < —0 and at x > +0 respectively, Ui are displacement of 
the neighboring tubes. (See Fig. Wp)- 

The equations for the neighboring tubes, coupled with tube "0" , are 

cm = K5{x)[{U+ - U{) + (t/_ - f/,)]. (3) 

One can see from the symmetry of equations that all Ui are the same. In what follows I 
put Ui = Ui. The solution can be found in a form 

f/_ = e'*=" + De"'^" (x < 0), f/+ = Ce'*^" (x > 0), 

U{ = Ae-''''' (x < 0), U^ = Ae^^^ (x > 0). (4) 



The time dependent factor exp{—iujt) is omitted here and below. Since the forces acting on 
the tubes are locahzed near the openings and the openings are supposed to be less than the 
wavelength, the displacements have discontinuities at point x = 0. The boundary conditions 
can be obtained by integrating Eqs.(PI21) over x. One gets that at x = 



dx dx 



The solution has a form 



\D 
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= r([/+ + t/_-2t/i), (5) 

f)TJ 

--^ = ^mr-U^)- (6) 



' (36i^2 + P)(49ir2 + A;2) 



(49^2 + p)(;^2 + ^2/36) 



(4227^2 + 36P) ' 



Due to the energy conservation 



12^2 + + = 1. (7) 

In the case of weak interaction k/K ^ oo one gets 

1^12 = 1,1^12 = 17112 = 0, (8) 

which means that the wave completely reflects from the terminational point. In the opposite 
case of a strong interaction k/K —>■ 0, 

|D|2 = 36/49, |C|2 = |A|2 = 1/49. (9) 

It will be important below that in the case of the strong interaction the structure of the 
operator C is irrelevant since I ignore the spatial derivatives in the boundary conditions 
Eq. ©. Therefore it does not matter whether or not this wave is acoustic and what spectrum 
it has. In fact, the numbers in Eq. @ is determined only by the geometry of the problem. 
For the geometry of Fig. ^ (c), the equations of waves are 

C{U) = K6{x){U2 + - 2U), (10) 
C{U^) = K6{x)iU-U^), (11) 
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where U is displacement of tube 0, C/^ and U2 are displacements of tube 2 at x > and 
X < respectively. The solutions have a form 

^- ^gite^^^g-ifex^ U+^Cxe'^'', (12) 

U^^Aie-'^'', U+^Aie^^''. (13) 

The boundary conditions at x = are 

dU~ dU'^ 

K{U^ + Ut-2Ul (14) 



dx dx 

dx 



Then 

Due to the energy conservation 2|Aip + + |Cip = 1. For the weak interaction K <^k 
one gets 

|Ci|' = 1,|^P = |L>i|' = 0. (17) 

In this case the wave is transmitted without any scattering. In the opposite case K ^ k 
one gets 

|Di|2 = |^|2 = |Ci|2 = l/4. (18) 

As in the previous case the result for the strong interaction is independent of the form of 

the operator C 

THERMAL CONDUCTIVITY 

Let us now calculate the TC of the perfectly aligned nanotubes. In the approximation of 
elastic scattering the heat flux Q along each row of the ahgned nanotube conserves because 
the waves generated in neighboring nanotubes due to scattering have zero total momentum. 
This leads to a conservation of Q along the row because in the theory of phonon thermal 
conductivity any relaxation of Q is the result of momentum loss. It is important, however, 
that at the points of scattering of both types, considered above, the numbers of phonons 
in each mode changes. Therefore the symmetric parts of the distribution functions in these 
points can be considered as in equilibrium with different temperatures for each point. Finally, 
I assume that between the scattering points of both types the propagation is ballistic. 
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FIG. 2: Part of the row of nanotubes with two cuts. Dotted hnes in the nanotubes correspond to 
the cuts of the neighboring nanotubes. The regions between them are considered as bahistic. 



In average every section of a nanotube can be divided into seven ballistic regions such 
that each boundary of the region corresponds to a cut in one of the six neighboring rows of 
the nanotubes. The part of one row is shown in Fig. |21 Since the energy flux is the same 
along the row, but the scattering is different the temperature intervals between neighboring 
boundaries are also different. To calculate the TC, I find the total temperature difference 
through all the nanotube at a given flux Q. 

Consider a region i of one nanotube and assume that each end of a region perfectly 
matches a thermal bath. The temperature difference of the left and right boundaries of the 
region is tj. Thermal flux produced in this region is Q = G£r)ti, where the function G{T) 
is called thermal conductance. It can be written in a form|l8l| 

where z = fuj/kBT, and the sum is over all monotonously increasing segments of spectrum 
^s{k), Zai and being the lower and upper boundaries of such segments. Here ks and h 
are the Boltzmann and the Plank constants respectively. To calculate the above integrals 
one should know the vibration spectra of nanotubes. They have been calculated previously 



17, hj, ab initio 



within different frameworks such as an empirical force constant model 
studies {20!. and tight-binding molecular dynamicsj^^. I use here the function G{T) calcu- 
lated by Yu. Gartstein[2^. His results for phonon spectra does not differ much from the 
previous ones though he used some original method. 

Suppose that the temperature decreases from left to right and the flux oriented in the 
same direction. The flux coming from the section to section 1 is QjCp, where |Cp is 
the transmission coefficient of the first scattering problem. Ballistic flux generated in the 
region 1 is —Gti, where G is ballistic conductance given by Eq. (jl9|) and ti is the negative 
temperature difference between boundaries 0-1 and 1-2 in Fig. 2. To get the total flux in the 
region 1 one should take into account reflection at the boundary 1-2. Simple calculations 
show that multiple reflections inside the same region do not change the result substantially. 



Finally, in the region 1 

Q = {Q\C\^-Gh){l-\D,\'). (20) 
Thus I can find the change of the temperature ti in the first region 



In the same way I find 



-«o = gf^47Tra-|C,n. (23) 



while 

^ G yi- \D\ 

The total change of the temperature AT through the segment can be calculated as 
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AT = J2u = -QB/GiT), (24) 



where 



B = { + ^^TT^ - 6|Cip - . (25) 



,l-|L>i|2 1-\D\ 

Since I suppose that the interaction is strong, coefficients in B are independent of the 
frequency. As a result I can express the TC k, through G{T). 

Flux per m^ is Qt = QN, where is the number of tubes in the bundle per m^. Assuming 
triangular lattice it is easy to get A^ = 4 x 10^''m~^ for (10,10) nanotubes 13, ^|. Using 



equation AT = LdT/dx, definition of TC Qt = —ndT/dx, and Eq.(j2H) one gets 

K = NLG{T)/B. (26) 

One can see from Eq.(jH))and Eq. (jl7|) that in the case of the weak interaction B ^ oo 
and K — > 0. That is because the jump of the wave through the cut is a bottle neck of the 
problem. Thus, assuming that interaction is strong, I get a maximum estimate for the TC. 

From Eqs. ()9ll8p one can find that 1/B = 0.0976, and Eq. ()26|) leads to the final result 

K = 0m7QG{T)LN. (27) 



where L is the length of one nanotube. 

Fig. 01 shows the results of Eq. (f?7|) at L = 1, 0.87, 0.7yU together with the results by Hone 
at a/. 4]. One can see that the theory refiects well enough both the magnitude and the 
temperature behavior. In fact, the only parameter here is the average length of a nanotube. 
The deviation at high temperatures is probably related to Umklapp processes. 
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FIG. 3: Thermal conductivity as calculated from Eq. ()27() for the (10,10) tube with L = 0.7 (dashed 
line), 0.87 (solid) and 1.0 (dotted) micron. The experimental data of Ref. J are shown by diamonds. 

CONCLUSIONS 

Finally, I present the maximum estimate of the TC of perfectly aligned nanotubes taking 
into account the scattering of phonons by the terminational points of the nanotubes. This 
estimate gives a quantitatively correct description of the thermal conductivity of aligned 
nanotubes as obtained experimentally by Hone et al. ^ | assuming that the length of segments 
is of the order of 1 /im. It follows from my results that the way to make thermal conductivity 
of the aligned nanotubes at room temperature larger than about 300 W/m K is to increase 
their lengths. Of course the TC will not increase indefinitely with L, as it follows from 
Eq. ()27j) . because sooner or later the mean free path due to other scattering processes will 
be smaller than L/7. However, some additional gain may be achieved with increasing L. 

I am grateful to V. Agranovich and A. Zakhidov for fruitful discussions. The work has 
been funded by DARPA subcontract SC0302 
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